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ABSTRACT 

The Phoenix Deep Survey is a multi- wavelength survey based on deep 1.4 GHz radio 
imaging, reaching well into the sub-100 /iJy level. One of the aims of this survey is 
to characterize the sub-mJy radio population, exploring its nature and evolution. In 
this paper we present the catalog and results of the spectroscopic observations aimed at 
characterizing the optically "bright" {R < 21.5 mag) counterparts of faint radio sources. 
Out of 371 sources with redshift determination, 21% have absorption lines only, 11% 
show AGN signatures, 32% are star-forming galaxies, 34% show narrow emission lines 
that do not allow detailed spectral classification (due to poor signal-to-noise ratio and/or 
lack of diagnostic emission lines) and the remaining 2% are identified with stars. For 
the star-forming galaxies with a Balmer decrement measurement we find a median 
extinction of Afja = 1.9 mag, higher than that of optically selected samples. This is 
a result of the radio selection, which is not biased against dusty systems. Using the 
available spectroscopic information, we estimate the radio luminosity function of star- 
forming galaxies in two independent redshift bins at 2 « 0.1 and 0.3 respectively. We 
find direct evidence for strong luminosity evolution of these systems consistent with 
LiAGHz oc {l + zf-'^. 
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Subject headings: galaxies: evolution — galaxies: starburst — radio continuum: galax- 
ies — surveys 

1. Introduction 

Deep radio surveys reveal, below a few milliJansky, a population of sources (the sub-mJy radio 
population) that is significantly different from the "classical" active galactic nucleus (AGN) powered 
radio galaxies typical at higher flux levels. This is clearly revealed in the normalized differential 
radio source counts, which shows a flattening at ~1 mJy that is not reproduceable by evolutionary 
models of the mJy and Jy radio population (Danese et al. 1985, 1987; Dunlop Sz Peacock 1990). 

Initial follow-up observations identified many of the sub-mJy sources with blue, often mor- 
phologically disturbed, galaxies (Kron, Koo, &; Windhorst 1985; Windhorst et al. 1985; Thuan & 
Condon 1987; Windhorst, Mathis Sz Neuschaefer 1990). Optical spectroscopy has revealed a large 
fraction of galaxies with spectra indicative of star formation (~ 1 — 100 Mq yr~^), similar to nearby 
galaxies detected by the InfraRcd Astronomical Satellite (IRAS), and extending to z ~ 1 (Frances- 
chini et al. 1988; Benn et al. 1993; Hammer et al. 1995; Georgakakis et al. 1999; Gruppioni, Mignoli, 
&; Zamorani 1999; Prandoni et al. 2001). Satisfactory models of the faint radio source counts re- 
quire the inclusion of an evolving population of such star-forming galaxies (Rowan-Robinson et al. 
1993; Hopkins et al. 1998, 1999) in addition to the AGN population. 

Interest in the sub-mJy population has increased significantly in recent years. Firstly, deep 
radio surveys are an efficient tool for compiling star- forming galaxy samples to z 1. This is 
essential to constrain the evolution of these systems, and to estimate the global star formation rate 
density out to cosmologically interesting redshifts. Unlike the optical and UV, radio wavelengths 
are insensitive to dust obscuration. Therefore, radio selection not only provides samples that are 
not biased against dusty starbursts, but also does not require the large and uncertain corrections 
for dust extinction that have plagued optical/UV studies (e.g., Afonso et al. 2003). Secondly, there 
is some overlap between the still mysterious sub-millimeter galaxy population, possibly dominated 
by dusty star-forming galaxies at high redshifts, and very faint (Si.4ghz ^ 200 /xJy) radio galaxies 
(Barger, Cowie, & Richards 2000; Carilli et al. 2001). The radio emission of these sub- mm sources 
is often the best available clue in refining the source position to allow optical identification (Blain 
1999). It also provides an initial, if coarse, redshift estimate (Carilli &; Yun 1999), which may be 
the only such estimate, since optical and near-infrared fluxes are often too faint for detection. 

Despite the potential of the faint radio population for cosmological studies, there is still limited 
information on the nature and evolution of these systems, primarily due to incomplete redshift 
information and small sample sizes. Haarsma et al. (2000) combined different deep radio surveys to 
constrain the cosmic star formation history to z ~ 2. Their star formation rate density estimates are 
in good agreement with IR-selected samples (Rowan-Robinson et al. 1997) and much higher than 
those obtained from UV and optical studies. These differences most likely arise from dust induced 
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biases affecting tfie latter samples. However, the combined sample of Haarsma et al. (2000) suffers 
from small number statistics (total of 77 sources) and inhomogenous selection criteria due to the 
different frequencies, radio flux density limits and spectroscopic identification rates of the individual 
samples used. A systematic study of a well defined sample of radio sources at microJansky levels 
is necessary to securely establish the nature and evolution of these systems. 

The Phoenix Deep Survey (PDS) is currently one of the largest area radio surveys reaching 
fiJy levels, comprising over 2000 radio sources within an area of over 4.5 deg^ to the flux density 
limit of about 60/xJy at 1.4 GHz (Hopkins et al. 1998, 1999, 2003a). Over the last few years, we 
have pursued a multiwavelength program to increase our understanding of the nature of the sub- 
mJy and ^uJy radio populations, and the evolution of star- forming galaxies (e.g., Georgakakis et 
al. 1999; Mobasher et al. 1999; Hopkins et al. 2000; Georgakakis et al. 2000; Afonso et al. 2001; 
Georgakakis et al. 2003; Afonso et al. 2003). In this paper we present the follow-up spectroscopic 
observations of the optically brighter {R < 21.5 mag) radio sources in this field, resulting in a large 
spectroscopic catalog of 371 redshifts that can be used for evolutionary studies. Spectroscopy of 
fainter optical counterparts {R > 22 mag) and photometric studies of yet fainter sources will be 
presented in forthcoming papers. 

The present paper is organised as follows: in Section 2 we describe the radio and optical 
observations, data processing and source classification using the spectroscopic data. Section 3 
presents the spectroscopic catalog, while in Section 4 the nature of the sub-mJy radio population 
and the obscuration of radio-selected star-forming galaxies are investigated. Section 5 explores the 
evolution of star-forming galaxies directly from the spectroscopic sample, through the determination 
of the radio luminosity function in two independent redshift bins at z pa 0.1 and 0.3. In Section 6 
the conclusions are summarised. 

Throughout this paper, and unless otherwise noted, we adopt a Hq = 70 kms~^Mpc~^, 
= 0.3, and J^a = 0.7 cosmology. 

2. OBSERVATIONS 

2.1. The Phoenix Deep Field radio observations 

The Phoenix Deep Survey (PDS) includes a deep 1.4 GHz survey made using the Australia 
Telescope Compact Array (ATCA) covering an area of 4.56 square degrees (the Phoenix Deep Field, 
PDF), selected to lie in a region of low optical obscuration and devoid of bright radio sources. A 
noise level of 12 /xJy RMS is measured in the most sensitive regions, using a resolution of ~ 6" x 12". 
Source detection resulted in a catalog of 2148 sources over the whole field, with integrated fluxes 
ranging from 59 ^Jy to 114mJy. Full details of the radio observations, data reduction and source 
catalogs are presented in Hopkins et al. (2003a). 
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2.2. Optical identifications 

Optical CCD imaging was carried out at the Anglo- Australian Telescope (AAT) in the R- and 
y-bands, as detailed in Georgakakis et al. (1999). The i2-band observations cover around 70% 
of the PDF, while the F-band covers a smaller region. These observations are relatively shallow, 
being complete to i? ~ 22.5 mag (Georgakakis et al. 1999). Deep wide field multicolour imaging 
(UBVRI) of the central part of the PDF has recently been acquired, and is presented in Sullivan 
et al. (2004). In the current paper, only the i?-band observations of these two surveys are used, to 
identify optical counterparts of faint radio sources. 

As described in Sullivan et al. (2004), a nearest-neighbour criterion was chosen to identify 
the best optical counterpart candidates. Using the i?-band optical images, the 4" radius region 
surrounding each radio detection was searched for an optical counterpart. If several possibilities 
exist in this region, the nearest source is assumed to be the counterpart, and a flag is set to indicate 
the existence of multiple candidates. To give some a posteriori quantification of the likelihood of 
any given optical identification, the reliability (TZ) was calculated for each identification, following 
the likelihood ratio method of Sutherland &: Saunders (1992). Optical identifications were found 
for over half (1181 out of 2148) of the sources in the radio catalog of the PDS. The different 
optical catalogs used prevent a simple conclusion regarding the optical identification rate being 
drawn. However, considering only the shallow i?-band observations, 46% of the radio sources are 
identified to R = 22.5 mag (the respective completeness level). When considering the deepest R- 
band pointing of Sullivan et al. (2004), optical identifications are found for 73% of the sources, 
down to R = 24.0 mag. These values are comparable to those found in previous work (Windhorst, 
Kron, k Koo 1984; Windhorst et al. 1985, 1995; Gruppioni, Mignoh, & Zamorani 1999). 



2.3. Spectroscopic observations 

Optical counterparts of radio sources in the PDS were observed spectroscopically in a series of 
observing campaigns. Initial data were obtained using multi-fibre spectroscopy at the Two Degree 
Field Spectrograph (2dF) at the Anglo-Australian Telescope (AAT) (Dec. 1996 and Sept. 1997) 
and slit spectroscopy at the ESO 3.6-m telescope (Oct. 1996). These observations were described 
in detail in Georgakakis et al. (1999). 

Further spectroscopy was obtained at the AAT with 2dF in September 1998, August 1999, 
October 2000 and July 2001. The 2dF is composed of two spectrographs, each with 200 fibers 
having a 2" diameter on the sky. Sources can be targetted over a two degree diameter region. This 
allows the simultaneous acquisition of up to 400 low to medium resolution spectra in a wide field, 
and is ideal for covering large survey areas such as the PDS. 

As in the initial 2dF observations, sources were separated into groups according to optical 
magnitude. The total exposure time varied from 1.5 h for the brightest souces to 3 h for the optically 
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fainter ones (i? ~ 21 — 22 mag), each comprised of separate half- hour exposures. The September 
1998 run was affected by poor weather, with only one half-hour exposure on the brightest objects 
being successful. Most of these sources were re-observed in later runs. 

Low resolution gratings were selected, being the most appropriate for efficient redshift iden- 
tification and classification. The 270R and 300B gratings were used on the 1998 run, while the 
270R and 316R gratings were selected for the 1999 and 2000 runs. The final 2001 run made use 
of the 300B grating only. The throughput of the 270R and 316R gratings produces an effective 
spectral window spanning the wavelength range 5000 < A < 8500 A, while that for the 300B is 
4000 < A < 7000 A. The spectral resolution is A/AA ~ 600 with a pixel scale of ~ 4.5 A pixel"^ In 
each run, objects were distributed over the available gratings to maximise the number of allocated 
fibers. 



2.3.1. Data processing 

The spectroscopic observations were processed using the 2dF data reduction (2dfdr,) package. 
Averaging of the reduced spectra for multiply observed individual objects was performed using the 
FIGARO package in the STARLINK environment. 

To estimate emission line fiuxes, a flux calibration was done in the following way. A relative 
flux calibration was performed using a response function provided by the 2dF team, calculated by 
estimating the 2dF efficiency in the V- and i?-bands and then interpolating using a second 
order polynomial. In order to have an estimate of the line-fluxes, the spectra were then calibrated 
using the measured total i?-band magnitude for the observed object. This provides an "aperture 
correction" which is usually small since most of the galaxies observed are spatially not much bigger 
than the fiber. For the more extended galaxies, however, the spatial undersampling can give rise to 
a bias in terms of classification (the fiber will be preferentially sampling the central regions of the 
galaxy) and in the line luminosity. Although the number of such objects is small in this sample, 
and does not affect the conclusions of the statistical analysis we present, any studies of individual 
bright objects in the PDS will need to take this effect into consideration. 

Redshifts were determined by visual inspection of the final spectra, through the identification 

of emission and/or absorption features. A quality flag, Q, was associated with each spectrum. A 
value of "A" denotes a high quality spectrum, with a firmly established redshift. A value of "B" 
indicates a low S/N spectrum, where a redshift was established using at least 2 spectral features. 
Spectra where only one spectral feature is identified, were assigned a Q value of "C" . 

Line parameters (flux and equivalent widths) were determined by Gaussian fitting to identified 
emission and absorption features using the splot package in iraf. Several factors affect the mea- 
sured line parameters. The low S/N of some of the spectra and the presence of strong atmospheric 
emission above 8000 A, not completely removed by the sky subtraction, affects the measurement 
of Ha at 2; > 0.25. Superposed stellar absorption, especially for the Balmer lines, will also affect 
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the measurements of emission in these lines. Stellar absorption at the wavelength of H/3 can be 
significant compared to typical emission line fluxes. This can be a major factor in biasing estimates 
of the true H/3 emission, and of derived quantities such as the Balnicr decrement for estimating 
obscuration. The average H/3 stellar absorption is around 2 A in star-forming galaxies (Tresse et 
al. 1996; Georgakakis et al. 1999; Hopkins et al. 2003b). Correction for this effect requires the 
knowledge of the H/3 absorption in each galaxy individually, something which can only be achieved 
with high-resolution and high S/N spectra. ConscqTicntly only the measured fluxes, including any 
underlying absorption, are presented here. When dealing with large samples, however, a statistical 
correction can be applied, as done below. 

In total, quality spectra were obtained for 371 optical counterparts of radio sources. Twenty- 
one sources were observed both in the earlier campaigns (Georgakakis et al. 1999) and the later ones, 
with seventeen of these having an Ha measurement. Redshift estimates agree within the measure- 
ment errors, while estimated Ha luminosities agree within 30%, the expected accuracy of the flux 
calibration of fiber spectra. Two galaxies, however, display a much larger difference (a factor of 2) 
in the Ha luminosity determination between the two runs. This discrepancy was traced to different 
fiber positions in extended (> 10") i bright galaxies, an effect previously mentioned. The positions 
observed were offset by 2", (comparable to the fiber size) between different spectroscopic runs, 
thus sampling different physical locations of the galaxy. As mentioned above, line measurements 
for these objects will not be representative of the whole galaxy. In the final catalog, measurements 
for objects with several independent determinations (such as these ones) were averaged together. 

Figure 1 shows the radio flux density and i?-band magnitude distributions for the 371 optical 
counterparts of faint radio sources with quality spectra. While the radio flux density range is well 
sampled, only for i? < 20 mag the optical magnitude range is signiflcantly covered. 

2.3.2. Spectral Classification 

The objects with spectral line measurements were classified according to their spectra. They 
were divided into the following categories: (1) absorption line systems, if showing only absorption 
features (Balmer lines, H+K 3950 A, 4000 A break, MgB 5175 A, NaD 5892 A); (2) Seyfert 1 ff pre- 
senting a broad line component in their emission line spectrum; and narrow emission line systems, 
either (3) Seyfert 2 or (4) star-forming, according to the properties of the emission lines detected. 
This separation of narrow emission line objects was done using the diagnostic emission line ratios 
[Oiii]A5007/H/3, [Sii]AA6716,6731/Ha, [Oi]A6300/Ha and [Nii]A6583/Ha, which are believed to 
identify the ionizing source responsible for the emission lines, ie., young hot (OB) stars or an active 
galactic nucleus (AGN) (Baldwin et al. 1981; Veilleux & Osterbrock 1987; Kewley et al. 2001). As 
the emission line pairs are close in wavelength, these ratios are insensitive to relative effects from 
dust obscuration and uncertainties in the absolute flux calibration of the spectra. At redshifts 
z > 0.3, however, the Ha line moves out of the spectral window of the observations, and the ratio 
[Oll]A3727/H/? must then be used (Rola, Terlevich & Terlevich 1997). While still sensitive to the 
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nature of the ionizing source, this ratio is more affected by such wavelength dependent effects. 

The classification assigned to a given source was that consistent with the majority of the line 
ratios available for that source. To minimize incorrect classifications due to the effect of dust 
obscuration, higher-z sources, with only [Oii]A3727/H/? and [Oni]A5007/H/3 line ratios measured, 
were considered star-forming only if they lie in the proper region of the diagnostic diagram after a 
correction for an obscuration corresponding to a Balmer decrement of Ha/H/3 = 6.0 (the average 
value for lower-2 sources in the present sample, as shown below) , using the Galactic extinction curve 
of Cardelli, Clayton k, Mathis (1989) with Ry = 3.1 - otherwise, the source remains unclassified. 

3. The spectroscopic catalog 

The Phoenix Spectroscopic Catalog was constructed using the full set of spectroscopic obser- 
vations of PDS radio sources described above. For higher redshift {z > 0.3) galaxies, where the Ha 
line could not be measured, an Ha luminosity was inferred using the [On]A3727 flux and an assumed 
ratio between [Oii]A3727 and Ha. This relation is known to be strongly luminosity and metallicity 
dependent (e.g., Jansen, Franx & Fabricant 2001), and this dependence has been suggested to be 
primarily an obscuration effect (Aragon-Salamanca et al. 2002). The relation employed here is the 
one found by Hopkins et al. (2003b) for radio-detected star-forming galaxies in the first data release 
of the Sloan Digital Sky Survey (SDSS, Abazajian et al. 2003), Fjon] = 0.23 x -Fka, which appears 
to be suitable for the present sample, as discussed in Afonso et al. (2003). 

The final spectroscopic catalog is presented in Table 1, and contains the following information: 

(1) radio source name; 

(2) Si. 4 GHz integrated flux density and its error; 

(3) right ascension (J2000) of the optical counterpart; 

(4) declination (J2000) of the optical counterpart; 

(5) i?-band magnitude (Vega-based system) of the optical counterpart; 

(6) reliability {TZ) of the optical identification, following Sutherland & Saunders (1992); 

(7) the measured [OII]A3727 equivalent width in A; 

(8) the measured H/3 equivalent width in A; 

(9) the measured Ha equivalent width in A; 

(10) observed Balmer decrement (Ha/H/?, uncorrected for any stellar absorption); 

(11) redshift z; 

(12) spectrum quality parameter, Q; 
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(13) spectral classification; 

(14) log of Ha (observed) luminosity in W; 

(15) log of 1.4 GHz luminosity in WHz~^. 

Out of the 371 optical counterparts of faint radio sources with a redshift determination, 78 
(21%) present only absorption lines in their spectra. Another 42 (11%) show signs of AGN activity 
(17 Seyfert 1 and 25 Seyfert 2) while 117 (32%) were classified as star-forming. For 126 (34%) 
sources a classification was not possible, due to the detection of too few emission lines and/or the 
low S/N of the lines detected. Eight sources (2%) were identified as stars. These can be cither true 
Galactic radio stars or foreground objects coincidcntally close to the radio position. Helfand et 
al. (1999) performed a search for radio stars using the FIRST^ radio survey, finding an extremely 
low detection rate at high Galactic latitudes (less than one star with radio emission per 200 deg^, 
for V < 10 mag). It is likely that the present sample, much deeper at both radio and optical 
wavelengths, may be detecting a higher density of radio-emitting stars. The reliability (TZ) of these 
identifications is in fact very high for 5 of the 8 sources, which suggests that most of these are 
actually Galactic radio stars. Very little is known on the luminosity function of radio stars and, 
hence, it is diflicult to estimate what number of such objects are expected in the high galactic 
latitude Phoenix Deep Field (see also the discussions in Kron, Koo, & Windhorst 1985; Oort & 
Windhorst 1985; Benn et al. 1993). All these sources are included in our spectroscopic catalog, but 
candidate stars are omitted from the analysis of the extragalactic radio population in the remainder 
of this paper. 



4. RESULTS AND DISCUSSION 

4.1. Nature of the sub-mJy radio population 

The present sample is used here to investigate the nature of the sub-mJy radio population for 
galaxies with optical counterparts. As shown in Figure 1, this will be representative of the R<20 
mag range. This includes a large number of sources (138) with radio flux densities below 200 /xJy, 
the spectroscopic properties of which are also explored. 

Figure 2 shows the radio flux density and i?-band magnitude distributions for the extragalactic 
radio sources in the present sample according to spectral type. The corresponding redshift distri- 
bution is shown in Figure 3. It is clear that the spectroscopic survey is particularly sensitive to 
the low and intermediate redshift (z < 0.5) objects. This will be to a significant extent due to 
the magnitude cut-off imposed by the spectroscopy. This can be seen from Figure 4, which shows 
the Hubble i?-band diagram for the present sample. Also shown are the predictions for M* early- 
and late-type galaxies (using the values from Nakamura et al. 2003), adopting the evolutionary 



^Faint Images of the Radio Sky at Twenty-cm (Becker, White & Helfand 1995). 
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scenario of Pozzetti, Bruzual, & Zamorani (1996) which make use of the population synthesis code 
of Bruzual & Chariot (1993). Extrapolating the broad relationship followed by faint radio sources, 
we can see that incompleteness due to the magnitude limit of the spectroscopic survey becomes 
important above z ~ 0.5. 

The increase of the star- forming galaxy population for sub-mJy fluxes (e.g., Benn et al. 1993; 
Georgakakis et al. 1999) is evident. While 11% of the sources with radio fluxes Si.4ghz > ImJy 
show a spectrum indicative of active star formation, at sub-mJy levels this number increases to 
nearly 40%. Table 2 summarizes the observed behavior, from the milliJansky down to the mi- 
croJansky level. The star-forming population in the present optical magnitude range appears to 
dominate at microJansky fluxes, while the number of absorption line and Seyfert galaxies drop 
steadily towards fainter radio fluxes. Unclassified objects, however, increase significantly at the 
faintest radio flux density levels. This can be attributed, at least partly, to fainter average optical 
magnitudes with decreasing radio flux density. As Figure 5 shows, the optical magnitude distribu- 
tion of the juJy radio sources appears to be lacking the brighter optical counterparts seen in the 
0.2 mJy < Si.4GHz < 1-0 mJy flux density range. Consequently, the median i?-band magnitude gets 
fainter by 0.7 magnitudes between the two radio flux density bins. Interestingly, this disapearance 
of bright optical counterparts for Si.4gHz < 0.2 mJy is reflected in a corresponding lack of lower 
redshift objects, when sub-mJy sources with fluxes above and below 0.2 mJy are compared (Fig- 
ure 5). This behavior could explain the higher fraction of spectra with lower S/N (fainter optical 
magnitudes) and where a classification is more difficult (sources at higher redshifts) . Deeper optical 
spectroscopy, currently underway, is thus necessary to perform a detailed study of the faintest radio 
sources. 



4.2. Obscuration of star-forming galaxies 

The present dataset can be used to investigate the obscuration affecting radio-selected star- 
forming galaxies. Previous studies (Hopkins et al. 2003b; Afonso et al. 2003, the latter using a 
preliminary sub-sample of the present catalog) have shown that radio selection, avoiding dust- 
induced biases, can give quite different results to optically-selected surveys. In particular, Hopkins 
et al. (2003b) derives a value of Ana = 1.2 mag for optically selected star-forming galaxies in the 
the first release of the Sloan Digital Sky Survey data (Abazajian et al. 2003), while a higher value of 
A-Ha = 1-6 mag is obtained for the subsample of galaxies with radio detection in the FIRST survey. 

Figure 6 shows the distribution of Balmer decrements (Ha/H/?) and derived emission-line 
obscuration values for Ha, Ahq, for star- forming galaxies in the PDS, using the standard Galactic 
extinction curve of Cardelli, Clayton & Mathis (1989) with Ry =3.1. Stellar absorption of the 
Balmer lines was corrected for assuming an average value of 2 A for the equivalent width (EW) 
of the H/3 absorption in star-forming galaxies (Tresse et al. 1996; Georgakakis et al. 1999), with 
a similar value (2.1 A) being used for the Ha line (equation (2) of Miller & Owen 2002). The 
corresponding distributions for the subsample of galaxies with radio detection in the FIRST survey 
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of Hopkins et al. (2003b) is also shown in this Figure by a sohd histogram, scaled down to allow 
the comparison of the shapes of the distributions. A median value of Aho = 1.9 mag is found, 
somewhat higher than the Ana = 1.6 mag estimated by Hopkins et al. (2003b). Although the 
spread in Ahq is large, and it might be tempting to think that there is little significant difference in 
these values, the distributions look different, the present sample being more uniformly distributed 
to higher values of Ahq. In fact, a KS test between the present sample and the radio-detected 
SDSS star-forming galaxies in Hopkins et al. (2003b) shows that the two arc not drawn from a 
single parent distribution at the 99.9% confidence level. The observed difference is likely to be due 
to the radio selection of the PDS, unbiased against highly obscured galaxies. We note, however, 
that the requirement of an optical spectroscopic measurement with Ha and H/3 may still result in a 
bias against the most obscured star-forming galaxies. It is likely that the true average obscuration 
value for star-forming galaxies is higher still. 



4.3. Nature of the unclassified radio population 

The existence of the unclassified radio population prevents a complete quantitative analysis 
of the distribution of sub-mJy and //Jy sources (for R < 20 mag) in star- forming and AGN sys- 
tems. For all of these sources, though, at least one emission line was observed. How can their 
characterization be taken further? 

A first approach is to study their radio luminosity. Even the most intense starbursts will not 
give rise to the highest radio powers observed in some AGN (-C/i.4GHz ^ 10^^^^^ WHz^-*^), and 
most star-forming galaxies will display radio luminosities substantially lower. The distribution 
of radio luminosity according to spectroscopic classification for the present sample is presented 
in Figure 7. The sources with Si.4ghz < 200 /iJy are represented by a dark cross-hatched his- 
togram. A weak bimodality exists between star-forming galaxies and absorption line systems or 
Seyfert galaxies. Star-forming galaxies cover the lO^''"^^'® WHz~^ range, with a median value 
of -^^i.4GHz ~ 10^^'^WHz^^. Sub-mJy Seyfert and absorption line galaxies cover a wider range 
^]^q20-26.5 -y^ jj^^i for the former, 10^*^"^^ WHz^^ for the latter), displaying a median value of 
-^i.4GHz ~ lO^^-^WHz"^. The faintest (Si.4ghz < 200 /iJy) radio population shows approximately 
the same behavior, although a lack of the highest luminosity objects is apparent. Also noteworthy 
is the existence of a number of absorption line and Seyfert systems with relatively low radio power 

(i^l.4GHz<1022WHz-l). 

With the lack of the star-forming systems above lO^^'^WHz"^, it might seem reasonable to 
classify otherwise unclassified systems with higher radio luminosities as AGN dominated systems. 
We note that this radio luminosity translates to a star formation rate of 1750 MQyr^^ (using the 
calibration from Bell 2003), and few galaxies are expected to have higher star formation rates. 
For the majority of the unclassified systems, though, the lack of a clear bimodality in the lumi- 
nosity distribution between star-forming and AGN systems makes a distinction impossible. Radio 
luminosity alone thus seems of very limited use for classifying sources in deep radio surveys. 
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The radio-to-optical flux ratio can also be investigated as a possible type discriminant for the 

unclassified objects. Using a large sample of milliJansky radio sources classified through their mid 
and far-infrared IRAS fluxes, radio morphology and radio polarimctry, Machalski & Condon (1999) 
showed that this parameter separates between normal/star-forming galaxies and AGNs (their Figure 
5). Following their procedure, we define the radio to optical flux ratio, ri.4, as 

n.4 = log(Si.4GHz//R) (1) 

where /r is the flux density at 6700 A calculated from the measured i?-band magnitude: 

fn = 3.0 X 106-'^-4*«(mJy). (2) 

The distribution of radio-optical flux ratio for the present sample, as a function of source type, is 
shown in Figure 8. 

The star-forming population shows a well defined distribution, very similar to the one found 
by Machalski &; Condon (1999). The absorption line and Seyfert galaxy populations, however, 
while showing a distribution with higher ri,4, are not as distinct from the star-forming population 
as observed by Machalski &: Condon (1999). The ri.4 distribution for AGNs in the present sample 
seems to be shifted to much lower values than those sampled in the Machalski &; Condon (1999) 
sample. For Si.4ghz < 200 /xJy this separation becomes indistinguishable, although here we are 
dealing with a relativelly small sample. 

Given the fainter radio and optical detection limits of the current survey, it is necessary to 
confirm that similar ranges of ri.4 are indeed being probed in each case. The sample from Machal- 
ski Sz Condon (1999) uses radio sources from the NVSS^ radio survey, with radio flux densities 
Si.4GHz > 2.5 mJy and optical counterparts brighter than i? ~ 18 mag. For the PDS, the limits 
are Si.4gHz > 0.059 mJy and R < 22 mag. The range of ri.4 spanned by both samples is shown in 
Figure 9. While there is a factor of ~ 50 separating the radio flux density limits sampled, there is 
an almost equivalent improvement in the optical magnitude Hmit for the PDS, allowing a similar 
range of ri.4 to be probed. Consequently, the bimodal distribution in ri.4 between AGNs and star- 
forming galaxies, if existent, should be seen in the present work. It is worth noting, however, that 
for Si.4GHz<200/xJy one would need to reach R ~ 26 mag to sample the highest values of ri.4. 

Being a flux ratio, ri.4 is sensitive to redshift effects. However, the variation of this parameter 
for the same galaxy within the redshift range of interest (z < 0.5, both in the present sample and 
in Machalski & Condon 1999) is, at most, a few tenths, due to the different rest-frame emission 
sampled. Furthermore, the redshift range sampled by Machalski & Condon (1999) and in the 
present sample is similar {z < 0.4 for the former, z < 0.5 for the latter). Considering a correction 



^NRAO VLA Sky Survey (Condon et al. 1998). 
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for the slightly higher redshift range sampled or, more directly, restricting the PDS sample to 
z = 0.4, does not change Figure 8 significantly. 

The differences in the radio-optical flux ratio distributions seem instead to be due to the 
different classification methods employed. The classification of Machalski & Condon (1999) use 
radio morphology (by looking for the existence of radio lobes or jets) and polarimetry to identify 
AGN galaxies, which guarantees a selection of high radio luminosity AGN with large values of 
ri.4, but misclassifying many low-power AGN (as indeed observed by Sadler et al. 2002). The 
spectroscopic-based classification for Phoenix sources, on the other hand, reveals low luminosity 
AGN equally well, and these show much lower values of ri.4. 

Figure 8 reveals some differences between the ri.4 distributions for star-forming and AGN 
galaxies, although they are not clear enough to easily disentangle the nature of the unclassified 
sources in this survey. The median ri.4 value for the unidentified sample is closer to the one found 
for AGNs than for star-forming systems, which could favor an AGN interpretation for the nature 
of the majority of these systems. There is another factor, though, that needs to be taken into 
account. As indicated by Barger, Cowie, & Richards (2000), there is some overlap between sub- 
mJy radio sources and galaxies revealed by mm and sub-mm surveys. These are thought to be 
dusty star-forming galaxies at higher redshifts z~l — 3, which would be revealed as sub-mJy and 
microJansky radio sources. As shown above, heavily obscured galaxies do exist in the PDS (see also 
Afonso et al. 2001). Radiation at optical wavelengths would be suppressed by dust, which would 
result in high radio-optical flux ratios, easily occupying any adopted "AGN regime". Figure 10 
shows some indication for this effect. Star-forming galaxies in the present sample are represented, 
with those that have a (Balmer decrement) obscuration measurement identifled. Among these, 
the ones with the highest obscuration (Aho > 3.0 mag) show a tendency for higher ri.4 values, as 
expected. Unfortunately, the number of galaxies with available Balmer decrement measurements 
is still too small to draw any reliable inferences. Larger surveys, such as the SDSS, would be ideal 
for investigating the relation between ri.4 and obscuration values. We note that the star-forming 
galaxies with the highest values of ri.4 in Figure 10 are at redshifts too high {z > 0.25) for a 
Balmer decrement measurement with the present set of observations. These could indeed be more 
obscured galaxies, since there is a selection of more luminous radio systems at higher redshifts, 
simply due to the flux limited nature of the present survey. As seen in Afonso et al. (2003) and 
Hopkins et al. (2003b), more luminous radio star-forming galaxies are also expected to display, on 
average, a higher obscuration. 



5. Evolution of star- forming galaxies 

Modelling of radio source counts at sub-mJy levels requires a population of evolving star- 
forming galaxies (Rowan- Robinson et al. 1993; Hopkins et al. 1998). Though attempts have been 
made to identify this evolution directly, through spectroscopy of optical counterparts (eg Mobasher 
et al. 1999), the number of radio sources with spectroscopic information has been too low to 
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result in direct evolution detection. Adding further information from photometric redshifts and 
classifications has been used to provide a successful indication of the evolution present in the sub- 
mJy radio population (Haarsma et al. 2000). 

The present sample, being more complete and homogeneous than previous investigations, al- 
lows for a direct test of evolution in the faint radio star-forming population. This is done here by 
calculating the radio luminosity function for star-forming galaxies in two independent redshift bins: 
"locally", defined here as z < 0.2, and at higher redshifts, 0.2 < ^; < 0.5. 



5.1. Radio Luminosity Function of star- forming galaxies 

To calculate the radio luminosity function for star-forming galaxies, only sources with this 
spectroscopic classification and radio flux densities in the range 150/LtJy< Si.4ghz < 3mJy inside 
the optically observed region were considered. The lower limit to radio flux prevents the faintest 
sources from being included in the luminosity function determination, since those sources are only 
detected in a small fraction of the survey area. Furthermore, at the faintest radio flux levels 
the unclassified sources become as abundant as the star-forming population, which could cause a 
significant error in the luminosity function determination. The i?-band magnitudes were limited to 
the completeness level of the optical survey, R < 22.5 mag. Incompleteness is then present due to: 
(1) only a fraction of radio sources having optical identifications and (2) only a fraction of these 
being observed spectroscopically and resulting in a redshift determination. The simplest correction 
for these biases can be done by defining completeness functions ?7i(Si.4gHz) and r]2{R), respectively 
the ratio of the number of radio sources with optical identification to the total number of radio 
sources and the ratio of the number of galaxies with a measured redshift to the total number of 
sources with optical counterparts in the sample. Figure 11 presents r7i(Si.4GHz) and r]2{R), showing 
the ~ 50% optical identification rate and the good spectroscopic coverage for R < 20 mag, where 
around 70% of the sources with optical identification have a redshift determination. 

The radio luminosity function was thus estimated for the star-forming sample described above 

as 

n ^ 

$(Ll.4GHz)=5] ^7 (3) 

1=1 ' 

with an associated error of 




with the summation over the number of galaxies in radio luminosity bins. Vmax is calculated, for 
each source, by 
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-fo Ai (1 + 2)2 VS!«(1 + 2P + S!a 

where is the luminosity distance, and the integral limits are zl = m8LK{zradio, ^optical, ^min) and 
z2 — rniii(^,.a(ijo' "^opticai' The quantities z^adio ZopUcal are, respectively, the redshifts 
corresponding to the bright radio flux (3.0 mJy) and optical r-band magnitude (R = 14 mag) limits 
of the sample, while z'^^^^^^ and z'^^-^i are the maximum redshifts at which the given source could 
still be detected (i.e., with Si.4ghz = 0.15 mJy and R = 22.5 mag). The quantities Zmin and z^ax 
are, respectively, the minimum and maximum redshifts covered by the galaxies in the radio sample 
considered. is the solid angle over which a given radio source could have been detected, given 
the noise level of the radio map, while ?7i(Si.4gHz) and r]2{R) are the completeness functions at 
the radio 1.4 GHz flux density and optical i?-band magnitude of that source (given in Figure 11). 
i?-band ii'-corrections were estimated using model spectral energy distributions for intermediate- 
type spirals, while at 1.4 GHz a spectral index of a = 0.8 (Si.4ghz oc z^~°), typical of star-forming 
galaxies, was assumed. 

Radio luminosity functions were estimated for "local" galaxies, with z < 0.2, and higher 
redshift galaxies, at 0.2 < z < 0.5. Table 3 lists the estimated luminosity function values in these 
redshift intervals. 

Figure 12 compares the derived luminosity functions with the local radio luminosity function 
for star- forming galaxies of Sadler et al. (2002) , in its parametric fit form (solid line) . While a good 
agreement is seen for the low redshift determination, the luminosity function for higher redshift 
galaxies is measurably higher. This is in agreement with the proposed evolution of the radio 
luminosity function (Rowan- Robinson et al. 1993; Hopkins et al. 1998; Haarsma et al. 2000). In 
particular, assuming a simple luminosity evolution for the Sadler et al. (2002) luminosity function 
in the form Li.4ghz oc (1 -|- z)'^, with Q = 2.7 ± 0.6 (Hopkins 2004) seems appropriate to model 
the observed difference (shaded region). In terms of star formation rate this corresponds to an 
increase from the local value of 0.018 Mq yr^^ Mpc^^ to 0.035 Mq yr~^ Mpc~^ at a median redshift 
of 0.32 (using the calibration from Bell 2003). This is also consistent with the evolutionary rates 
derived in previous work (Hopkins et al. 1999; Haarsma et al. 2000; Seymour et al. 2004), but arises 
directly from a pure spectroscopic sample of radio selected star-forming galaxies. A good coverage 
of the sub-mJy radio population to an optical magnitude of i? < 21 mag seems to be necessary 
to start directly detecting evolution in the faint radio star-forming galaxy population. Extending 
the optical magnitude range covered by the spectroscopy will result in strong constraints on the 
evolutionary status of the sub-mJy radio star-forming population, a work currently underway in 
the Phoenix Deep Survey. 
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6. Conclusions 

A catalog of spectroscopic observations of faint radio sources in the Phoenix Deep Survey 
has been presented and used to investigate the nature of the sub-mJy radio sources. Down to 
R < 21 mag we find a substantial fraction of star- forming galaxies (36%), with a smaller fraction 
of (presumably) AGN host galaxies (27%), and a large fraction of emission line galaxies for which 
no classification was possible (35%). Commonly used criteria (-Li.4gHz and radio-optical flux ratio 
ri.4) have been unable to establish the nature of the unclassified population, essentially due to the 
sensitivity of this survey to low-power AGNs. 

The 1.4 GHz luminosity function for star-forming galaxies in the present sample, estimated 
using the spectroscopic information directly available, shows good agreement with previous deter- 
minations of the radio luminosity function for z < 0.2. At higher redhifts (0.2 < z < 0.5), the 
estimated luminosity function shows an evolution over the local determination, compatible with a 
simple model of luminosity evolution of Li.4ghz cx: (1 -|- z)'^, with Q = 2.7 it 0.6. 
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Table 1. The PDS spectroscopic catalog 
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Table 1 — Continued 
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Table 2. Galaxy populations at faint 1.4 GHz flux levels 



Galaxy type 


Sl.4GHz > 1-0 
(mJy) 


0.2<Si.4GHz<1.0 

(mJy) 


Sl.4GHz <0.2 
(mJy) 


star-forming 


11% 


37% 


36% 


absorption 


47% 


20% 


12% 


Seyfert 


14% 


14% 


7% 


unclassified 


28% 


29% 


45% 
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Table 3. Radio (1.4 GHz) Luminosity Function for star-forming galaxies 



log-Ll.4GHz 


(1)2=0-0.2/ r \ 


ff)2=0.2-0.5/' r \ 


(W Hz-1) 


(Mpc~^ dex-^) 


(Mpc"3 dex-^) 
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23.1 




(1.2 ±0.6) X 10"^ 


23.5 




(1.3 ± 0.7) X 10-5 
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Fig. 1. — Radio (1.4 GHz) flux density distribution (a) and i?-band distribution (b) for the 371 
optical counterparts of faint radio sources with a spectroscopic redshift determination. The open 
histogram in (b) represents the i?-band magnitude distribution for the full sample of optical iden- 
tifications of the Phoenix Deep Survey. 





Fig. 2. — Radio (1.4 GHz) flux density distribution (a) and i?-band distribution (b) for the extra- 
galactic sources in the present survey according to spectral type. 
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Fig. 3. — Redshift distribution for the faint radio galaxies in the sample according to spectral type. 
Eleven Seyfert 1 and one unclassified galaxy with 1.0 < z < 3.1 are not shown in this Figure. 
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Fig. 4. — /2-band magnitude vs. redshift for the sources shown in Figure 3, according to spectral 
type. The sohd and dashed Hnes represent, respectively, the predicted variation of aparent mag- 
nitudes for M* early- and late-type galaxies, according to the evolutionary scenario of Pozzetti, 
Bruzual, & Zamorani (1996). The dotted line shows the R = 21.5 mag level, indicative of the 
spectroscopic target selection magnitude limit. 
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Fig. 5. — Magnitude (left) and redshift (right) distribution for different types of faint radio sources 
with redshift determination in selected ranges of 1.4 GHz flux density. The median magnitude and 
redshift for each radio flux interval is indicated. 
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Ha/H/3 A„„ (mag) 

Fig. 6. — Distribution of Balmer decrements (corrected for stellar absorption) of star- forming galax- 
ies in the present sample (left) and derived (emission-line) obscuration values at the Ha wavelength 
for the same galaxies (right) . The solid histograms represent the (scaled down) distributions found 
by Hopkins et al. (2003b) for radio-detected SDSS star-forming galaxies. See text for details. 
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Fig. 7. — Radio luminosity distribution of sub-mJy radio sources according to their spectroscopic 
classification. The median radio luminosity value is indicated for each class. Sources with 1.4 GHz 
flux densities below 200 /xJy are also represented (dark cross-hatched histogram). 
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Fig. 8. — Radio-optical flux ratio distribution of sub-mJy radio sources, according to their spec- 
troscopic classification. The median radio-optical flux ratio value is indicated for each class. The 
dark cross-hatched histogram represents sources with 1.4 GHz flux densities below 200 /xJy. 
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Fig. 9. — Radio-optical flux ratio parameter space covered by the study of Machalski &; Condon 
(1999, open squares) and the sub-mJy sources described here (filled triangles). Diagonal lines show 
the locus of different values for the ri.4 parameter, as indicated, while the dashed horizontal line 
shows the radio 1.4 GHz flux density limit of the Phoenix sample. 
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Fig. 10. — Distribution of radio-optical flux ratio for star-forming galaxies. Galaxies with a (Balmer 
decrement) obscuration measurement are identified (hatched histogram) and, among these, those 
with the highest obscuration measurements, Ahq > 3 mag (solid histogram). 
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Fig. 11. — Completeness functions used to correct for incompleteness in the determination of the 
1.4 GHz luminosity function for star-forming galaxies in the PDS. r7i(Si.4GHz) is the ratio of the 
number of radio sources with optical identification to the total number of radio sources within 
the survey area considered. r]2{R) represents the ratio of the number of galaxies with a measured 
redshift to the total number of sources with optical counterparts in the same sample. Each identified 
star- forming galaxy contributes to the 1.4 GHz luminosity function correcting for the incompleteness 
given by ?7i(Si.4ghz) and ri2{R) at that galaxy's radio flux density and i?-band magnitude. 
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Fig. 12. — Radio (1.4 GHz) luminosity function for star-forming galaxies in the PDS, locally {z < 
0.2, dots) and at intermediate redshifts (0.2 < z < 0.5, triangles). The parametric fit form for the 
luminosity function derived by Sadler et al. (2002) is also shown, both in its derived form for local 
galaxies (solid line) and after applying a simple luminosity evolution in the form Li.4ghz oc (l + z)'^, 
with Q = 2.7 ± 0.6 as derived by Hopkins (2004) (shaded region). 



